Vortex configurations of bosons in an optical lattice 
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The single vortex problem in a strongly correlated bosonic system is investigated self-consistently 
within the mean-field theory of the Bose-Hubbard model. Near the superfluid-Mott transition, 
the vortex core has a tendency toward the Mott-insulating phase, with the core particle density 
approaching the nearest commensurate value. If the nearest neighbor repulsion exists, the charge 
density wave order may develop locally in the core. The evolution of the vortex configuration from 
the strong to weak coupling regions is studied. This phenomenon can be observed in systems of 
rotating ultra-cold atoms in optical lattices and Josephson junction arrays. 
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Vortices in charged or neutral superfluids (SF) can be 
created by applying an external magnetic field |l| or by 
rotating the system Q. Properties of vortices are essen- 
tial for understanding superfluidity Q. Generally, vor- 
tices are topological defects in the SF order parameter 
with a 2n phase winding around the core. The SF or- 
der is suppressed near the core in an area characterized 
by the size of the healing length £. As a result, other 
orders competing with the SF order may develop there. 
In the context of high T c superconductors, it was first 
proposed in the SO (5) theory that the vortex core is 
antiferromagnetic in underdoped samples 4]. This the- 
oretical prediction has now been confirmed by numerous 
experiments 5] . Using the physics of the vortex core to 
study competing orders of doped Mott insulators is an 
important new direction in condensed matter physics 0. 

Quantized vortices have been observed in rotating sys- 
tems of dilute ultra-cold alkali atom gases 0, , drawing 
intense attentions both experimentally and theoretically 
HHIHEEElHIlil- Most discussions up to now 
have been based on the Gross-Pitaevskii-Bogoliubov (G- 
P-B) equations, which assume that the particle density 
is given by the square of the amplitude of the SF order 
parameter. Thus, the minimum particle density is al- 
ways located at the core, shown experimentally as a dark 
region in vortex imaging. 

Near the SF-Mott insulator (MI) transitions, the G-P- 
B method ceases to work well. Theoretically, many inves- 
tigations based on the Bose-Hubbard model are available 
now 0, in the absence of rotation. The MI phases 
occur at the commensurate (integer) fillings when the 
hopping amplitude t is small compared with the Hub- 
bard repulsion U . Transitions into the SF phase are 
triggered by either increasing t or changing the chem- 
ical potential. In the presence of the nearest neighbor 
repulsion W, the insulating CDW phases appear at half- 
integer fillings. However, the vortex configuration in this 
region has not been fully investigated. Experimentally, 
tremendous progress has been made in realizing the SF- 
MI transition in ultra-cold atomic systems in optical lat- 
tices 0, 113, as proposed in Ref. plj. It would be 
interesting to further generate vortices in such systems. 



Alternatively, Josephson junction arrays and granular su- 
perconductors in magnetic fields are also possible systems 
to study such vortices, where Cooper pairs behave like 
composite bosons. 

In this article, we address the vortex configuration near 
the SF-MI and SF-CDW transitions. In the former case, 
the vortex core is close to the MI phase with a nearly 
commensurate filling while the SF dominates in the bulk 
area. In the latter case, the superfluid vortex with the 
CDW core behaves as a "meron" topological defect of the 
3- vector psuedospin order parameter. These strong cou- 
pling configurations evolve to the weak coupling G-P-B 
vortex smoothly as t/U increases. Although the discus- 
sion below is for the rotating neutral bosonic system, 
it is also valid for the charged system in the magnetic 
field. Theoretical predictions obtained in this paper can 
be tested in these systems. 

We study the 2D Bose-Hubbard model extended by a 
nearest neighbor repulsion W in the rotating frame 
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with /i the chemical potential, m the boson mass, Q = Clz 
the rotation angular velocity, the position r related to 
the axis of rotation, the vector potential A due to the 
Coriolis force, V c f the centrifugal potential, and V ex the 
trapping potential. U is scaled to 1. For the charged 
bosonic system in the magnetic field B, A = (e* /2fic)B x 
r instead and the V c f term is absent. 

This model can be solved self-consistently by the mean- 
field (MF) approximation 0, |2l|, by decoupling the 
inter-site terms as 
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FIG. 1: (a) The ground state phase diagram for Bose- 
Hubbard model at W = and Q — 0. Lines of equal particle 
densities are plotted with (N) = 2.5 ~ 1.5 from top to bot- 
tom, (b) The SF order parameter (a) versus (N) at different 
values of t/U = 0.1 ~ 0.01 from top to bottom. 




FIG. 2: Vortex configurations near the SF-MI phase bound- 
ary at t/U = 0.02 with bulk particle density (N) = 1.95 (in 
(a) and (b) ) and (N) = 2.05 (in (c) and (d)). The left hand 
side (a) and (c) are the particle density distributions; the right 
hand side (b) and (d) are the SF order distributions. 

n(fi)n(fj) « n(f f i )(n(f f j)) + (i j) - (n(r*t)}(n(rj-))- (2) 

(a(fj)), (n(fj)) are expectation values on the MF ground 
state, which is approximated by a product wave func- 
tion of the form \^)g = n^ 4 \ip(fi)). We cutoff each 
single-site Hilbert space up to 10 particles which is suf- 
ficient for experimental values (N) w 1 ~ 3 19]. The 
SF and CDW order parameters are defined as (a(fj)) 
and (— ) Vi ((n(fi)) — (N)) respectively, where (N) is the 
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FIG. 3: The evolution of vortex configurations at (N) = 1.95 
as increasing t/U. SF amplitudes in (a) and particle densities 
in (b) are shown along a path cut from (10,20) to (30,20) in 
the 40 x 40 system. t/U = 0.02 ~ 0.08 from bottom to top in 
(a) and top to to bottom in (b). 

bulk average particle density. This simple MF theory de- 
scribes the SF-MI transitions and extrapolates well into 
the intermediate coupling region 0,0]. For example, 
the relation of the SF order (a) v.s. t/U from the MF 
theory is in good agreement with that of the Monte-Carlo 
simulations at both commensurate (integer) and incom- 
mensurate fillings [22]]. 

Before discussing the vortex problem, it is helpful to 
recall the well-known phase diagram when f2 = |l7j 
(Fig. ^ (a)). From the equal density contours we can 
see an approximate particle-hole (p-h) symmetry with 
respect to (N) = 2 near the phase boundary: fi drops 
or increases with increasing t/U when (N) > 2 or < 2 
respectively. These can be viewed as particle or hole SF, 
but this difference no longer exists as t/U becomes large, 
where fi drops with increasing t/U for both cases, i.e., 
the system evolves from the strong to weak coupling re- 
gion. In Fig. H](b), SF order vs. (N) is shown at different 
values of t/U. The SF order parameter increases mono- 
tonically with (N) at large t/U, while the commensurate 
filling suppresses the SF order prominently at small t/U, 
eventually leading to the MI phase. At t/U ps 0.08, the 
suppression disappears even at (N) — 1 , which marks the 
cross-over into the weak coupling region. 

We study the single vortex problem in a 40eto x 40an 
(an being the lattice constant) system around which the 
circulation of A is 2-7T and thus fl = h/ (2mL 2 )(L = 40an) 
correspondingly. The rotation center is located at the 
center of the central plaquette. We further simplify the 
problem by dropping the V ex and V c f terms, since they 
behave smoothly near the center of the trap and they can 
cancel each other if the is close to the trap frequency. 

Two typical vortex configurations near the SF-MI tran- 
sition are shown in Fig. [2] (a) ~ (d) with t/U = 0.02, 
where the particle densities in cores are maximal or min- 
imal respectively. The vortex core is located at the center 
of the plaquette with the reduced SF order on the sites 
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FIG. 4: The evolution of vortex particle density distribution 
at fixed t/U = 0.03 (a) and 0.02 (b) with varying bulk values 
(N) along the same path in Fig. £3 (a) From top to bottom, 
(N) = 2.3 ~ 1.7. (b) From top to bottom, (N) = 1.9 ~ l.f . 



nearby. Roughly speaking, the square of the superfluid 
amplitude, |(a)| 2 , is proportional to \(N) — Nq\ near the 
transition, where Nq = 2 here is the nearest commensu- 
rate value. The local particle density in the core should 
be closer to the commensurate value in order to suppress 
the SF order. As a result, the particle density reaches 
a maximum or minimum when the bulk density (N) is 
slightly smaller or larger than the commensurate value. 
The former case can also be understood as the vortex 
of the hole superfluid, where the hole density goes to a 
minimum at the core. This contrasts with the case of 
fcrmionic superfluidity, where the Cooper pairs are bro- 
ken into normal particles in the core with total particle 
density almost unchanged, and also with the case of the 
weak coupling bosonic system where the only possibil- 
ity is the depletion of the core particle density. As we 
approach the vortex core from outside, the hopping pro- 
cess is frustrated and thus effectively t/U becomes small 
because of the phase winding of the SF order. As a re- 
sult, the vortex core is driven closer to the MI state than 
the bulk area. We also check the vortex configuration 
with integer value of (N) = 2 at the same value of t/U, 
where the particle density distribution is uniform with 
suppressed SF order in the core. 

We further discuss the evolution of the vortex core con- 
figuration from the strong to the weak coupling regions. 
The bulk particle density is fixed at (N) = 1.95 with in- 
creasing t/U as shown in Fig. |2| The SF order increases 
and the healing length £ decreases away from the SF-MI 
boundary along this direction of evolution. On the other 
hand, the weak coupling expression £/ao = \/t/{U{N)) 
states that £ increases with increasing t/U . Thus we can 
infer that £'s minimum value appears in the intermedi- 
ate coupling region. The vortex core with extra particle 
density survives until t/U « 0.05 ~ 0.06, after which 
it crosses over gradually to that with depleted particle 
density at large t/U. This transition agrees with the be- 
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FIG. 5: (a) Phase diagram of the extended Bose-Hubbard 
model with W/U = 0.1 around (N) — 1.5. Lines of equal 
particle densities are plotted with (N) = 1.6 ~ 1.4 from top 
to bottom. The CDW (b) and SF (c) order distributions in the 
vortex configuration near the CDW phase with t/U = 0.023 
and bulk average (N) = 1.5. 



havior of the SF order vs. (N) when the system is not in 
rotation. In that case, at the same filling level, we check 
that the suppression due to the commensurate filling also 
vanishes around a similar value of t/U. When t/U is 
larger than 0.06, the vortex configuration is already sim- 
ilar with the weak coupling case. We also check the above 
evolution with fixed (N) = 2.05. The feature of SF or- 
der is similar with that in Fig. [3] (a) , and the minimum 
particle density is always located in the core. When t/U 
is less than an intermediated value around 0.05, the core 
particle density is close to the commensurate value 2.0. 
When t/U grows larger, it drops further. This also agrees 
with the evolution picture from strong to weak coupling 
physics. 

Another evolution from the particle-like vortex core to 
the hole-like core with fixed t/U in the strong coupling re- 
gion and varying (N) is shown in Fig. 01(a). We choose 
the region close to the MI phase with (N) = 2, where 
the approximate p-h symmetry is valid, with t/U fixed 
at 0.03 and (N) varying from 2.3 to 1.7. As a result, 
the difference between the core particle density and the 
bulk value changes from negative to positive as the bulk 
density passes (N) = 2. This point is the closest one to 
the tip of the MI phase along the evolution, where the 
minimum of the SF order and the maximum of £ also 
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lie. In Fig. 01(b), such evolution is shown along the path 
t/U = 0.02 connecting two neighboring MI phases with 
(N) = 2 and 1. The core configuration is close to the 
MI phase with (N) = 2 on the top, and becomes close 
to the MI phase with (N) = 1 at the bottom. Around 
(N) = 1.55, the density distribution becomes almost uni- 
form. This point is the maximum of the SF order and the 
minimum of £ because it is the farthest point from the 
MI, which is just opposite to Fig. E|(a). In both Fig. 0Ja) 
and (b), the tendencies to MI phase in the vortex core are 
strong when the average (TV) is close to an integer, and 
becomes weaker as (N) away from commensurate fillings. 

Next we turn on the nearest neighbor repulsion W/ U = 
0.1. In the absence of rotation, CDW phases appear be- 
tween two neighboring commensurate MI phases at half- 
integer fillings when t/U is small 23]. At large values of 
t/U, SF phases stabilize as usual. Between the CDW and 
the SF phases, the mean field theory gives a small area of 
the coexistence of CDW and SF order, i.e., the supersolid 
phase (2j| . These are shown in Fig. [5] (a) with the equal 
density lines around (N) =1.5. It is well known that 
the hard-core boson model can be mapped into spin 1/2 
XXZ model in a magnetic field, and thus the CDW and 
SF orders can be unified in a 3- vector pseudospin picture. 
With the releasing of the hard core constraint, the above 
mapping is still approximately valid in the sense that the 
spin up and down states correspond to two nearest inte- 
ger number states on each site. The vortex configuration 
at t/U = 0.023 with bulk particle density (N) = 1.5 is 
shown in Fig. [S] (b) and (c). The SF order dominates 
outside the core, while the CDW order develops together 
with the suppression of the SF order in the core. In the 
pseudo-spin picture, this is a kind of topological defect 
called "meron" , where the pseudo-spin pointing along z 
axis at the origin gradually changes to lying in the x-y 
plane with a winding number 1 around the origin when 
far from it. This vortex configuration also evolves to the 
weak coupling one at large value of t/U with the disap- 
pearance of the CDW order in the core. This situation 
is similar to the behavior of the superspin in the SO(5) 



theory of the antiferromagnetic vortex core[4], |5|. 

In previous experiments, the vortex core with size 
2£ w 0.4/xm is too small to observe directly by opti- 
cal methods. A time-of-flight expansion is needed be- 
fore optical absorption imaging |g. On the optical lat- 
tice, the vortex core is larger. The typical core size in 
our calculations is estimated at 5 ~ 6 lattice constant ao 
(a — 0.426 /jm in Ref. [l9|), i.e. about 2 /j,m. The reso- 
nant probe laser beam can be focused to this size at the 
level of current technology (24|. Thus without turning 
off the trap, it is possible to image the core particle den- 
sity distribution non-destructively by scanning the probe 
laser beam. It would be interesting to find the anoma- 
lous vortex with the maximum particle density in the 
core. Another possible realization is the Josephson junc- 
tion array. The non-uniform charge distributions in the 
vortex configuration result in electric fields [2^|. Thus 
it is possible to determine the filling in the vortex core 
with respect to the outside by measuring electric field 
distributions. 

In summary, we have studied vortex structures of the 
strong coupling boson systems. Near the SF-MI transi- 
tion, the vortex core is more strongly coupled compared 
to the bulk area and is thus closer to the MI phase with 
suppressed SF order. The particle density in the core can 
be either the maximum or the minimum of the whole sys- 
tem, always approaching the nearest commensurate den- 
sity of the Mott insulator. Near the SF-CDW transition, 
a superfluid meron-like vortex is found with a CDW core. 
All of these strong coupling vortex configurations evolve 
to the conventional weak coupling one as t/U increases. 
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